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We demonstrate that in two-dimensional noncentrosymmetric s-wave superconductors under ap- 
plied magnetic fields for a particular electron density, topological order emerges, and there exists a 
zero energy Majorana fermion mode in a vortex core, which obeys non-Abelian statistics, in anal- 
ogy with p x + ip y superconductors, the Moore-Read Pfaffian quantum Hall state, and the gapped 
non-Abelian spin liquid phase of the Kitaev model. 
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Recently, there is considerable interest in emergent 
topological states of many-body quantum systems char- 
acterized by a topologically nontrivial structure of the 
Hilbert space, i.e. "topological order" [lj]. In certain 
classes of topological states in 2+1 dimension, quasipar- 
ticles are non-Abelian anyons 0, S, H S, EH, 0, B, @, E3| • 
The essential feature of the non-Abelian statistics is that 
the exchange of particles is described by unitary opera- 
tions in the multi-dimensional space, which is spanned 
by the basis of the degenerate many-body ground state. 
Thus, the state depends on the order of the multiple ex- 
change processes of particles. The possible realization of 
non-Abelian statistics in real systems has been exten- 
sively studied so far in connection with the v = 5/2 
fractional quantum Hall (FQH) state, and the vortex 
state of chiral p x + ip^ superconductors (superfiuids) 
0, H i, H i, 0, El lH. These topological states 
are of interest also in the context of quantum comput- 
ing, since the non-Abelian anyon can be utilized as a 
decoherence-free qubit, and potentially applied to the 
construction of fault-tolerant topological quantum com- 
puters [H, El, 14, El. In this paper, we present an- 



other candidate of a topological phase allowing the ex- 
istence of non-Abelian anyons, which can be realized in 
strongly noncentrosymmetric (NC) s-wave superconduc- 
tors. This topological phase belongs to the same class 
as those of the Moore-Read (MR) Pfaffian FQH state 
0> Px + iPy superconductors and the gapped non- 
Abelian spin liquid phase of the Kitaev model ||, El El • 
In NC superconductors, the asymmetric spin-orbit (SO) 
interaction which breaks inversion symmetry plays im- 
portant roles in various exotic superconducting proper- 
ties [13, EES E|. in our proposal, the asymmet- 
ric SO interaction combined with an external magnetic 
field yields the nontrivial topological state for a particu- 
lar electron filling. 

We consider type II NC s-wave superconductors with 
the Rashba type SO interaction in two dimension. We ne- 
glect the parity-mixing of triplet components of Cooper 
pairs due to the asymmetric SO interaction [IT]], because 
the inclusion of this effect does not change the essential 



part of our argument. For a concreteness, we define our 
model on the square lattice, though the following con- 
sideration do not reply on the particular choice of the 
crystal structure. Then the model Hamiltonian is 

k.a k,aa' 

-J2lH^-u + h.c] (1) 

fe 

v=± fc fc 



Here, cl (c/. CT ) is a creation (an annihilation) opera- 
tor for an electron with momentum k, spin a. The en- 
ergy band dispersion is e k — —2t(cosk x + cosk y ) — /i. 
The second term of eq.fl]) is the Rashba SO interaction 
with jCo(fc) = (sin k y , — sin k x , 0). Eq.([2]) is expressed 
in terms of the chirality basis which diagonalise the SO 
term. The energy band is split into two parts by the 
SO interaction: e k ± — s k i a\Co(k)\. The gap func- 
tion in this basis is odd parity, and possesses the mo- 
mentum dependence A±(k) — Ar]^(k) with rj±{k) — 

— (Cox ± i£>Qy)/ J £q x + £-o y , which is, importantly, sim- 
ilar to that of p x + ip y superconductors. Thus, for 
A <C Ep, we can exploit the same argument as that 
applied to p-wave superconductors 2^, 23| , and find that 
in the mixed state with vortices parallel to the z-axis, 
there is a zero-energy quasiparticle state of a vortex core 
which is described by a Majorana fermion. However, in 
this case, the existence of Majorana fermions in vortices 
does not directly lead to the non-Abelian statistics of 
them, because there are two bands (fffe M , (J, = i), each 
of which contributes to a zero-energy mode with a dif- 
ferent band index. The existence of two different species 
of Majorana fermions in a single vortex implies that the 
sign change of the fermion operators under the braiding 
of two vortices, which is a hallmark of the non-Abelian 
statistics 0, Q (see below) , is canceled. To eliminate this 
unwanted multiplicity of Majorana fermions, we tune the 
chemical potential as /i = —At for which the Fermi level 
crosses the T point in the Brillouin Zone (BZ). In this 
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situation, there are still two bands near the Fermi level 
m the model ©: one from £fc_ with a finite Fermi mo- 
mentum, and the other in the vicinity of the T point 
which is given by the Dirac cone. To generate the mass 
gap in the Dirac cone, we introduce the Zeeman coupling 
HbH z Y,k( c k] c kl - c ii c *a)- Tne magnitude of the gap is 
of the order ^H z . Then, there is only a single energy 
band which crosses the Fermi level. Let us assume 
that H z is sufficiently smaller than the orbital depairing 
field -fforb • The Pauli depairing effect due to H z is negligi- 
ble for a 3> hbHz, A [19(. Under these circumstances, we 
can integrate out contributions from quasiparticles with 
a gap ~ /iB-fTz at the T point, and obtain the low-energy 
effective Hamiltonian TL c s for the single band supercon- 
ductor. Taking account of the fact that the Zeeman field 
H z induces the inter-band Cooper pairing between the 
(+)-band and the (— )-band [2l|, we obtain 
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]ek-al_a k - - ^[A_(fc)4_at k _ + h.c.]. (3) 



£fc- 



Here the renormalized energy band is £&_ 
£ with £ = H 2 A 2 /a 2 \C(k F )\ 2 m Q , \C{k)\ = 

^C 2 x + C 2 Qy + ^ B H 2 z /a 2 , and m = (4a 2 |£(fc F )| 4 + 

A 2 (Cl x + £ 2 ly ))/2a\C(k F )\ 3 . m a is an energy gap 
of quasiparticles in the (+)-band in the vicinity of 
the Fermi momentum of the (— )-band. The super- 
conducting gap function is A_(fc) = aAf]+{k) where 
a = 1 + H 2 A 2 /2a 3 \C{k F )\ 3 m 0l and T}+(k) = -(C 0x + 
iCoy)/\C{k)\. The above expression ([3]) is valid only 
in the vicinity of the Fermi momentum kp defined by 
£fc F _ = (not Sk F - =0). In many type II superconduc- 
tors, it is typical that for H z < H orh , [i-qH z < A. This 
implies that the gap generated by the magnetic field at 
the r point might be smaller than the superconducting 
gap, which may invalidate the approximation used in the 
derivation of j3j. However, our argument on the low- 
energy vortex core states which we are most concerned 
with is applied only to energy scale < A 2 /Ep. Therefore, 
the effective Hamiltonian ([3J is applicable for our pur- 
pose, as long as the condition A 2 /E F < hqH z < H OI b is 
satisfied, which can be fulfilled in ordinary experimental 
situations. 

The topological order of the model © 
clearly manifests in the Chern number which 
is, for the Hamiltonian of the form H c b — 
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- fc _)cr Al i; At (fc)(a fe _,aL fc _) t , defined 
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where E = (E x (k), E y (k), E z (k))/\E(k)\. The integral 
of ([4]) is taken over the whole BZ, while the expression of 
([3j is derived for k in the vicinity of k F . Nevertheless, we 
can consider the Chern number of the model ([3]) by re- 
interpreting eq.Q as a lattice regularized version of the 



low-energy effective theory, and extending the fc-space in 
which the model ([3J is defined to the entire BZ. Then, 
the numerical evaluation of J\f for the Hamiltonian ([3]) 
gives Af = 1. Therefore, the model (J3J is classified as the 
same topological class as those of the MR state, spinless 
Px + ipy superconductors, and the gapped non-Abelain 
phase of the Kitaev model. The existence of the Zeeman 
field H z in the model ((3j is important for this topological 
characterization, because it does not only break time- 
reversal symmetry, but also ensures the differentiability 
of E X:V (k) for Eq.([3|) which is singular at k — for H z = 
0. 

The Chern number M = 1 implies the existence of 
zero-energy Majorana fermion modes in vortices which 
obey the non-Abeian statistics, as in the case of p x + ip y 
superconductors d, H, To demonstrate this, we pro- 
ceed to solve the Bogoliubov de-Gennes (BdG) equations 
for the model ([3J with a single vortex inserted parallel 
to the z-axis. For simplicity, we switch to the contin- 
uum model replacing the energy band Sk of eq.{T]) with 
e' k = k /2m — fi and Co with C' = (k y , —k x , 0). Further- 
more, we assume that the gap amplitude A(r) vanishes 
inside of the vortex core, and is equal to a constant A 
outside of the core, and A <C E F . Then, in the vicinity 
of the Fermi surface, the BdG equations corresponding 
to the model ((3j with a single vortex are 
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where = (u(r),v(r)), P = -{d x +id v ), P f = -P*, 
and Ao = aA(r)/\C(k F )\. The BdG equations © are 
equivalent to those of spinless p x + ip y superconduc- 
tors except that there are the £o-terms in the diagonal 
components, which can be formally absorbed into the 
shift of the Fermi momentum kp — > kp — Eq/v f . Thus, 
the solution of (J3J) is given by * = e~ ieo/vi "^ p+ip with 
typ+ip the eigen function of the BdG equations for spin- 
less p x + ip y superconductors, and there exists a zero en- 
ergy mode inside the vortex core which is separated from 
the first excited state by a gap of energy size A 2 /E F 
[22I [23l | . The Bogoliubov quasiparticles for this zero 
energy state are described by a Majorana fermion field 
7 = J drlu^aLir) + v(r)a-(r)], since (v*(r), u*(r)) = 

(u(r),v(r)) for e = 0. Here a}_(r) = J^k a \- e ~ lkr ■ 

To confirm the above prediction, we apply numeri- 
cal analysis directly to the BdG equations for the tight- 
binding model |T]) without referring to the low-energy 
effective theory ([3J). The energy spectrum and the eigen 
functions of the BdG equations were calculated for the 
model (H|) with a vortex located at the center of the sys- 
tem on the square lattice with open boundaries. In this 
calculation, we assume that the GL parameter is so large 
that the Zeeman field H z is approximated to be uni- 
form, and the spatial dependence of the superconduct- 
ing gap function due to the vortex is taken into account 
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FIG. 1: The density of quasiparticles on the 37 x 37 xy-plane 
for e = 4.136 x W~ 4 t (top), e = 1.289 x 10~ 3 t (middle), and 
£ = 1.579 x 10~ 3 t (bottom). The left and right panels are, 
respectively, the plots of |M(r)| 2 and |i;(r)| 2 . 



only in its phase for simplicity. The topological prop- 
erties which we are concerned with are not sensitive to 
these approximations. We set parameters as [i = —4t, 
a = t, A = OMt, and ^ B H Z = 0.04. In FIGQ] shown 
are the spatial distributions of the density of Bogoliubov 
quasiparticles for several low-energy states calculated for 
the lattice size 37 x 37. The lowest energy state with 
e = 4.136 x 10~ 4 i is dominated by a vortex core state, 
which can not be the Caroli-de-Gennes-Matricon mode of 
the conventional s-wave superconductors, because for our 
choice of the parameters, the Fermi energy is Ep = 0.25£, 
and A 2 /Ep = O.Oli. Also the lowest energy level de- 
creases toward zero as the system size increases. Thus, 
we identify the lowest energy state with the zero energy 
mode. Furthermore, we find the low-energy edge states at 
the boundaries; e.g. for e = 1.289 x 10~ 3 , 1.579 x 10~ 3 . 
The edge state with energy < A is a concomitant of 
the zero energy vortex core state, which is in accordance 
with the Chern number Af = 1 @ . Taking these observa- 
tions into account, we can conclude that the NS s-wave 
superconductor with a magnetic field for the particular 
electron filling is in the topological state. 

As was proved by Ivanov, vortices with zero energy 
Majorana modes obey non-Abelian statistics Here 
we summarize some fundamental properties of the non- 
Abelian statistics relevant to the following argument 
[E @, M- Two Majorana fermion fields ji and jj in 
two vortices can be fused into a complex fermion field 
tp = (ji + i-fj)/2. The fermionic state described by ip 
is occupied or unoccupied. The fusion processes of vor- 
tices are analogous to the operator product expansions of 
the Ising conformal field theory with the central charge 
c = 1/2. The fusion rules for three particle states, a vor- 



tex (denoted as a), a fermion occupied state (tp), and 
an unoccupied vacuum state (1), are a x a = 1 + ip, 
ip X a = a, and ip X ip — 1. The degeneracy of the 
ground state with 2n vortices obtained from the fusion 
rules is 2 n ~ 1 . The braiding of vortices is given by a 
unitary operation in this multi-dimensional degenerate 
space. The braiding operator for vortices at i and j 
is R™ = CT exp(f 7j7l ) = ff ffltp(if (2^ty - 1)). Here 
a is a phase factor called a topological spin. Under 
the exchange of two vortices, the Majorana fermion op- 
erators are transformed as, ji — > R c " 7 ji(R a ' a ')^ = 7 j, 
7i — > R <T<T jj(R a ' 7 ) Jf — — 7j. The minus sign in the second 
transformation rule is essential for non-Abelian statistics. 
The topological spin 9 a in R aa is associated with the con- 
formal spin of the primary field for vortices which is the 
chiral Ising spin field, because of the bulk-edge correspon- 
dence of anyonic particles: 9 a — e l27T (h -h ) _ e %f _ j± s a 
result, when two vortices are fused into ip, the braiding 
of them yields the phase factor R™ = e 1 ? 1 , while when 
the fusion gives the topological charge 1, the phase factor 
due to the braiding is R\ a = e~ 1 ^ . The braiding rules 
for the other particle types are R°"v = —i and R^ = — 1. 

We, now, discuss the feasibility of the experimental de- 
tection of the non-Abelian statistics. One promising ap- 
proach is to use the two-point-contact interferometer pro- 
posed in the context of the FQH state pllll. 12 1. In the 
superconducting state considered here, this experiment 
is applicable only to the thermal transport. According 
[l2 |. the interference term of the edge heat 
depends on the parity of the total number of 
vortices n in the bulk. J lnt for odd n is much smaller than 
jmt £ or even n ^ though both of them do not exhibit the 
dependence on a magnetic flux $, because the cr-particle 
is neutral. This parity dependence characterizes the non- 
Abelian statistics. Another possible experiment is a bit 
indirect but simpler. It uses a disk-shaped system with 
which two heat baths are attached at the boundary, (see 
FIGl2Ja).) For this geometry, as in the case of p x + ip y 
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superconductors [24|, the energy spectrum of the edge 



state depends on the parity of the total number of vor- 
tices n in the bulk. For even n, the lowest energy state 
has a gap of the order A/kpL where L is the length of 
the boundary. Although the gap is small for a sufficiently 
large system size, it is nonzero, and thus the quasiparti- 
cle corresponding to this edge mode is a complex fermion 
interacting with This Bogoliubov quasiparticle is cat- 
egorized as the same particle type as the ^-fermion in the 
bulk, because, in the limit that the two vortices merge to- 
gether at a position r, the resulting ^-particle is nothing 
but the Bogoliubov quasiparticle with a nonzero energy 
Q. For odd n, the low-energy edge state is a Majo- 
rana fermion mode, and can be fused with an unpaired 
Majorana fermion in the bulk resulting in the ^-state 
or the 1-state. The phase accumulated by the current 
flow of the edge a- or V'-particles encircling the bulk n 
vortices is obtained from the square of the braiding op- 
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FIG. 2: (a) Experimental setup for the thermal transport 
measuremnt of the edge state. A magnetic flux $ penetrates 
in the bulk. Two heat baths are attached at the boundary so 
that the longer path between the two heat baths, C2, encircles 
almost all vortices in the bulk, (b) Superconductor-metal- 
insulator junction. 

erator (R ah ) 2 (a,b — <r,ip). When the temperature T\ of 
the heat bath 1 is smaller than the temperature X2 of the 
heat bath 2, the chiral edge heat current flows mainly in 
the path C2 (anticlockwise direction) encircling the bulk 
vortices. In this case, for even n, the edge heat current 
carried by the -0-particles exhibits a usual dependence on 
i.e. Jcvon ~ Em=i A m cos(27rme$//ic), while, for odd 
n, using the fusion rules and the braiding rules mentioned 
above, we obtain J^ d ~ Ylm=i ^ 4m C0S7Tm : where the 
i?4 m -term corresponds to a trajectory winding around 
the boundary loop Am times. Thus J"^ d is much sup- 
pressed. On the other hand, when T\ > T 2l the edge 
current flows mainly in the path Ci, less affected by <fr. 
In this case, the dependence of J mt on the parity of n is 
weaker than the case of T\ < T 2 . These observable ef- 
fects can be utililzed for the detection of the non-Abelian 
statistics. 

An advantage of NC s-wave superconductors over 
Px + ipy-wave superconductors and the v = 5/2 FQH 
state is that the gap energy scale of the former can be 
typically much larger than those of the latter. Note that 
the superconductivity in NC systems needs not to be a 
bulk phenomenon. Let us consider the junction between 
an s-wave superconductor and a metallic thin film placed 
on an insulating substrate, (see FIG. 2(b).) The thin film 
must be sufficiently clean so that the mean free path is 
larger than its thickness. In this system, inversion sym- 
metry is broken, and an asymmetric potential gradient 
perpendicular to the interface is introduced. We can use 
a material with a high transition temperature such as 
MgB 2 (T c ~ 39K) for the superconductor 25]. Then, the 
proximity effect induces s-wave superconductivity in the 
2D NC system realized in the thin film. If the Fermi en- 
ergy of the film Ep is much smaller than that of the bulk 
superconductor, the energy gap in the vortex core A 2 /Ep 
for the proximity-induced NC superconductor can be rel- 
atively large. The strength of the asymmetric SO inter- 
action can be controlled by changing the substrate or 
applying a perpendicular voltage on the film. Although 
electrons should experience strong SO scatterings at the 
interface, the transition temperature and the gap of the 



s-wave pairing state are not affected by them. Also, note 
that the Majorana fermions in vortices of the NC super- 
conductors do not require a half quantum vortex, i.e. a 
texture of the d-vector, because our system is essentially 
regarded as spinless. In this sense, zero energy Majorana 
states in NC s-wave superconductors are more realizable 
than in spinful p-wave superconductors. 

In conclusion, NC s-wave superconductors under mag- 
netic field have a topological order for a particular elec- 
tron filling, and can be playgrounds for the non-Abelian 
anyons. Although we consider only the Rashba SO in- 
teraction here, our argument can be easily generalized to 
other asymmetric SO interactions. 
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